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Abstract. We study a model of multiple-field DBI inflation leading to mixed form of pri- 
mordial non-Gaussianity, including equilateral and local bispectrum shapes. We present a 
general formalism based on the Hamilton-Jacobi approach, allowing us to go beyond slow-roll, 
combining the three-point function for the fields at Hubble-exit with the non-linear evolution 
of super-Hubble scales. We are able to obtain analytic results by taking a separable Ansatz 
for the Hubble rate. We find general expressions for both the equilateral and local type non- 
Gaussianity parameter /nl- The equilateral non-Gaussianity includes the usual enhancement 
for small sound speeds, but multiplied by an analytic factor which can lead to a suppression. 
We illustrate our results with two scenarios. In the first model, previously found to have 
detectable local non-Gaussianity, we find that the equilateral signal is not sufficiently sup- 
pressed to evade current observational bounds. In our second scenario we construct a model 
which exhibits both a detectable equilateral /nl and a negative local /nl- 
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1 Introduction 



Inflation is widely believed to be responsible for the specific set of initial conditions on which 
the hot big bang relies. Whilst a compelling mechanism, a consistent model proves elusive 
however (see [1-3] for reviews), which is in part due to the limited information available 
in the two-point statistics of the primordial density perturbations. Potential non-Gaussian 
signatures have therefore become an increasingly popular observable with which to discern 
between otherwise degenerate models [4], particularly since impending observations are set to 
improve by at least an order of magnitude [5] . It is therefore important to try to understand 
the correspondence between inflationary dynamics and the different forms of non-Gaussianity 
(see [6, 7] for recent reviews). 

Non-Gaussianity can be produced by inflation in a number of distinct ways. For ex- 
ample, by converting between entropy and adiabatic modes during 1 multiple-field inflation 
[15-17], the curvature perturbation £ can evolve on super-horizon scales [18, 19]. Such non- 
linearities can in principle produce local-type non-Gaussianity [20-42], where the effect is 
associated with a turn in the trajectory and is often enhanced by violations of slow roll. Al- 
ternatively, single field models with non-standard kinetic terms provide an alternative source 
of non-Gaussianity (see [43, 44] and references therein). Often motivated by string theory, 
the models we are concerned with have a characteristic sound speed c s , where c s = 1 in the 
canonical case. 2 As a result, equilateral non-Gaussianity can be produced by the interactions 
of quantum fields on sub-horizon scales. This is the case in Dirac-Born-Infeld (DBI) inflation 
[46, 47], in which a probe D-brane moving along the radial direction of a warped throat drives 
inflation. 

More generally however, it is expected that both equilateral and local contributions will 
be relevant in models characterised by non-standard kinetic terms and multiple-field dynamics 
[48-60]. For example, in previous work [59] we studied a multiple-DBI model, akin to that 
of [56-58], as a concrete example of multi-component inflation with non-standard kinetic 
terms. Using the 6N formalism, we tracked the super-horizon evolution of perturbations 
using the field fluctuations at horizon exit and the subsequent background trajectory. With 
the adoption of a sum separable Hubble parameter, as in [29], we were able to treat the 
two-field case both analytically and beyond slow variation to calculate the local contribution. 
Moreover, by considering inflation in the tip regions of two warped throats, we illustrated 
that rapidly varying sound speeds can produce large local type non-Gaussianity during a 
turn in the trajectory. 

Our previous work did not include the equilateral contribution produced on sub-horizon 
scales however, which is what we intend to address in this paper. Whilst this contribution is 
dominant in the single field case, the introduction of multiple fields can alter this conclusion 
through the conversion of entropy and adiabatic modes (see, for example, [48-50]). In this 
paper we again consider the multiple-DBI model as in [59] and compute the full third order 

1 Non-linearities can also develop through the conversion between entropy and adiabatic modes after infla- 
tion. For example, the curvaton mechanism [8-10] and modulated reheating [11, 12]. See [13, 14] for recent 
realisations of such scenarios in string theory. 

2 Although generally referred to as the sound speed this is technically the phase speed of fluctuations. See 
[45] for further clarification of this point. 
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action using the Arnowitt-Deser-Misner (ADM) formalism [61]. After considering the leading 
contributions in slow variation and small sound speeds, we calculate the three point function 
for the field fluctuations at horizon exit using the path integral approach [62]. Thereafter, we 
implement the 5N formalism and assume a separable Hubble parameter as in our previous 
work to calculate, fully analytically, the combined local and equilateral contributions to the 
bispectrum of the curvature perturbation, giving one of the few explicit examples of models 
characterised by both contributions (see [60, 63] for alternatives). Finally, as a first step 
towards assessing the viability of such a signal, we apply our results to two specific cases. 

The outline of the paper is as follows. We begin in section 2 by briefly reviewing the 
multiple-DBI model and introducing and reformulating the relevant quantities using the 6N 
formalism. Thereafter, we use the path integral method to calculate the three point function 
of field fluctuations at horizon crossing in section 3, having first derived the third order action 
for this scenario. We then use the SN formalism to present the corresponding equilateral 
non-linearity parameter. In section 4 we assume a separable Hubble parameter to combine 
this result with the local contribution found in our previous work, giving the total three point 
function for the curvature perturbation. We briefly assess the feasibility of such a signal by 
studying two specific examples in section 5. Finally, we conclude in section 6. 

Throughout this paper we use the (—,+,+, +) metric signature and set Mp = c = 1, 
where Mp = J— is the reduced Planck mass. Capital latin indices label scalar fields and 



any summation is explicit. Greek indices label space-time co-ordinates whilst lower case latin 
indices label spatial co-ordinates only, where the Einstein summation convention is adopted. 
Finally, commas denote partial derivatives and over-dots represent derivatives with respect 
to cosmic time. 

2 Multiple-DBI inflation and non-Gaussianity 

We begin this section by briefly reviewing the multiple-DBI model and the relevant obser- 
vational quantities, paying particular attention to non-Gaussianity. Thereafter, we use the 
5N formalism to re-write these expressions and investigate their evolution on super-horizon 
scales. This section is intended to be relatively brief, since full details can be found in [59]. 

2.1 Background evolution in multiple-DBI 

Multiple-DBI inflation is encompassed by the following action 



where Pj is a function of the single scalar field (pi and kinetic function Xj = —-^Q^^l^l^v, 
whilst the potential V is a function of the set of scalar fields (p = {(pi, <p%, 4> n }- We model 
inflation through n probe D3 branes descending n distinct warped throats glued to a compact 




(2.1) 
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Calabi-Yau manifold in type IIB string theory. 3 First considered by [56-58] to investigate 
the effect of multiple sound horizons on equilateral type non-Gaussianity, the corresponding 
expression for Pi is 



Pl = j m [l-}Jl-2fWX I ), (2.2) 

where parameterises the warped brane tension of throat I and is a function of <pj only. 
The corresponding homogeneous equations of motion are given by 



i, = -2c(')ff,, (2.3) 
Here c^ 1 ' is defined as the sound speed associated with the field I 



= = , 1 (2.5) 



and we have used (2.3) to eliminate Xj such that f^' remains a function of <pi whilst H, 
, V and <pj generally depend on the collection of fields eft. Notice that the above equations 
of motion are written in Hamilton- Jacobi form [65, 66] in which the Hubble parameter is 
written as a function of the scalar fields, taking precedence over the potential. This is more 
suited to the case of non-trivial sound speeds and will allow us to consider departures from 
slow variation. To define slow variation we introduce the following parameters 



/ i v 7 

s - Hcd)-/^ 3 ~ 2-j ,J c (j) fj ' W 



where we require e < 1 for inflation. Slow variation is defined as eW,r/^,s^ <C 1 and we 
shall state explicitly when this additional restriction is required. 

3 Note that this is distinct from the case of a single brane descending a warped throat along both radial 
and angular coordinates, as in the first example of multiple- field DBI [64]. 
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2.2 Perturbations and non-Gaussianity 



We characterise the scalar degree of freedom in the primordial density perturbations by 
introducing the primordial curvature perturbation on uniform density hypersurfaces x l ) 
(see [67, 68] for explicit definitions). The two and three-point correlation functions then 
define the power spectrum P^ and bispectrum respectively 



(C kl Ck 2 ) = (27T) 3 P C (fc 1 )5 3 (k 1 +k 2 ), 

(C kl Ck 2 Ck 3 ) = (27rfB ( (k 1 , k 2 , A; 3 )<5 3 (k! + k 2 + k 3 ), 



(2.7) 
(2.8) 



where £k is the Fourier transform of (, kj are comoving wavevectors and <5 3 is the three 
dimensional Dirac delta function. For a Gaussian £ the two-point function completely de- 
fines the statistics of the field. Signatures of non-Gaussianity are encoded in the connected 
contributions to higher order correlators, such as the three-point function. To parameterise 
the deviation from Gaussianity we introduce the fc-dependent non-linearity parameter /nl, 
given by the ratio of the bispectrum to a combination of power spectra 

6, , , v_ B c (k 1 ,k 2 ,k 3 ) 

g/NL(«i,A&,*3) - P dkl ) P< (k 2 ) + Pdh)Pdh) + Pd^PdksY { ] 

By assuming a scale invariant dimensionless power spectrum = Jr^P({k) the above can 
be written as 



pNL(ki,k 2 , k 3 ) - ^—^ ^-^2 • (2-10) 

We now adopt the 5N formalism [18, 70-73] to evolve £ on super-horizon scales using 
only the field fluctuations at horizon exit and the homogeneous field evolution thereafter. 
To facilitate this we make two restrictions on the background dynamics. First we demand 
that the sound speeds are comparable whilst observable scales exit during inflation, such 
that a- 1 ' ~ c* for all / during this interval. 5 Horizon exit 6 therefore equates to evaluating 
a quantity when c*k = a±H+. To simplify the spectrum of field fluctuations we also assume 
slow variation during horizon exit. 

Given these restrictions, we use the separate Universe approach [18, 71, 72, 74] to 
identify the curvature perturbation £ with the difference in the number of e-folds between the 
perturbed (N) and homogeneous background (Nq) universes, evaluated between an initially 

4 Here we adopt the sign conventions of [29] for ease of comparison. See [69] for a summary of the various 
conventions in the literature. 

5 This need not necessarily be the case, as in [57, 58] for example. 

6 We will use the more succinct term 'horizon' as opposed to 'sound-horizon' since there should not be any 
ambiguity as, in this context, the sound- horizon is the only relevant scale. 



-4- 



flat hypersurface (e.g. shortly after horizon exit) and a final uniform density hypersurface 
tf (e.g. early in the radiation dominated epoch). This allows us to calculate the relevant 
quantities (e.g. /nl) at time tf given the field fluctuations at time t* and the homogeneous 
field evolution between these times. For example, the dimensionless power spectrum can be 
written as [73, 75] 



where Nj is with respect to the field I at horizon exit. Here we have defined the dimensionless 
power spectrum of scalar field fluctuations at horizon exit using the two point function 

(<% kl <50 Jk2 ) = (2Tr) 3 5 IJ ^V*S 3 (k 1 +k 2 ), p+=(El) , (2.12) 



k\ * ' ' * \2it J 

where we use slow variation at horizon exit and 5 IJ is the kronecker delta symbol. Similarly 
the three-point function is given by 



UK 



\ NiNjN KL (5(f> Ikl 5(f>j k:i (6(t>K*5(f> L )k 3 ) + 2perms ), (2.13) 



2 

IJKL 



where in this case * denotes a convolution and 'perms' denotes cyclic permutations over the 
momenta. Neglecting the connected part of the four-point function and using Wick's theorem 
to rewrite the four-point functions as products of two-point functions, the latter term can be 
written as [73, 75] 



^ N ,i N ,J N ,KL{5^i\<i- s _8<l>j\<i 2 {5^K^8(j) L \ s ) +2perms = 

IJKL 

(2vr) 3 A^Vl E N jKjN,ij <5 3 (kx + k 2 + k 3 ), (2.14) 

such that the bispectrum becomes 

B c (h,k 2 ,k 3 ) = 4vr 4 P f 2 ^f (ji 3 l(h,k 2 ,k 3 ) + 5/<g) . (2.15) 
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Inspection of /nl(&i> &2, ^3) = /nl(^i> &2 5 ^3) + /nl shows that there are two distinct con- 
tributions to the bispectrum. Adopting the notation of [21], the fc-independent parameter 7 
/^L is due to non-linea 
contribution, given by 



is due to non-linear behaviour in £ on super-horizon scales and is referred to as the local 



6 ,(4) = EijNiNjNjj 
5 /NL 



/NL - /v- M 2 X2 • </- 10 J 



The result of our previous work was to provide an analytic expression for this parameter in 
the subset of cases described by a sum-separable Hubble parameter, in which the derivatives 
Nj and Njj can be fully evaluated. We neglected the contribution from the /c-dependent 

parameter f^{ki, k2, k^) however, which is due to the intrinsic non-Gaussianity of the dip 1 , 
produced by quantum field interactions on sub-horizon scales. 

The aim of this paper then is to explicitly calculate the equilateral contribution and so 
arrive at the total expression for the bispectrum. Inspection of (2.13) shows that this requires 
two distinct steps. In the following section we use the path integral method to first calculate 
the three point function of field fluctuations at horizon crossing. Thereafter we use the 8N 
formalism to find the equilateral non-linearity parameter of the curvature perturbation in 
this scenario. 



3 The equilateral contribution 

In this section we use a standard prescription to calculate the three point function of field 
fluctuations at horizon crossing, analogous to calculations in [49, 50, 62, 78]. We begin 
by presenting the third order action of field fluctuations for a more general scenario in the 
spatially flat gauge, before restricting ourselves to slow variation and small sound speeds 
around horizon exit in the multiple-DBI case. Thereafter, we use this result and the path 
integral formalism to find the three point function of field fluctuations at horizon exit and in 
turn an expression for fj?). 

3.1 The third order action 

To calculate the three point function of field fluctuations we begin with the general action 
(2.1). For a spatially flat Friedman-Robertson- Walker (FRW) Universe, the background 
equations of motion are given by 



7 For scale dependent cases see [76, 77]. 
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3H 2 = £ [tfPiM - Pi 
i 



+ V, 



(3.1) 
(3.2) 



which here we write in the conventional form, as opposed to the Hamilton- Jacobi form in 
section 2.1. The Klein-Gordon equation, which is not independent of (3.1) and (3.2), is given 
by 



Pu = mp I>Xl j>i + Pi,xr<j>i + Pi,xAi + Vj. (3.3) 

Progressing to perturbations about the homogeneous background, we construct the third 
order action by recasting (2.1) using the Arnowitt-Deser-Misner (ADM) formalism [61]. This 
will be useful since the lapse function N and shift vector A* become Lagrange multipliers 
under variation. This, along with an appropriate choice of gauge, will simplify the task 
of isolating the physical degrees of freedom when we consider perturbations. Until then, 
however, we stress that the equations remain exact with no choice of gauge. The ADM 
metric is given by 



ds 2 = -N 2 dt 2 + hij {dx i + PPdt) (dx j + N j dt) , (3.4) 

where hij is the spatial 3-metric. In terms of this metric, the action (2.1) and kinetic term 
become 



= - J d 4 xVh NR® + NKijK ij - NK 2 + 2N ^ Pj - 2NV 
^ I 



2N 2 



2 1 



(3.5) 
(3.6) 



where K = Kf, is the three dimensional Ricci scalar and indices are raised and lowered 
using the spatial metric. Kij is the extrinsic curvature, given by 



K„ 



2A7h 



\j + N j\i ~ h ij ) . 



(3.7) 



where u denotes the covariant derivative with respect to the spatial metric. To derive the 
energy and momentum constraint equations we vary the ADM action (3.5) with respect to 
the lapse function N and shift vector N l respectively 
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# (3) - K i:j K ij + K 2 - 2 



P I,X! 2 
N 2 1 



Pi 



+ V = 0, 



(3.8) 



(3.9) 



where for notational convenience we have introduced vj = $i—N % $i^. To solve the energy and 
momentum constraints we consider a first order 8 expansion of the inhomogeneous quantities 
about a spatially flat FRW background 



N 



91 = 91 + dpi, 
1 + a, N t -- 



((l-2V)%+2£|y), 



(3.10) 



where we consider scalar perturbations only and 5ij is the Kronecker delta symbol. This 
presents n + 4 scalar degrees of freedom: 6(f>i, a, (3, ip and E. We can eliminate two degrees 
of freedom by adopting the spatially flat gauge, whereby ip = and E = 0, such that 



91 = 91 + 091, 
N = l + a, Ni = P\ u 



h 



2 x 

a dij, 



(3.11) 



leaving n + 2 scalar degrees of freedom. Note that for notational convenience we drop the 
overbar on homogeneous quantities for the remainder of this section. To eliminate two 
further degrees of freedom we substitute the above into the constraint equations (3.8) and 
(3.9), giving algebraic equations for a and /3 

r 



a 



Pi,xM 

-2H—' 



(3.12) 



/3' 



I* 2H 



on ^ 



Pi,iXj<i> 2 i8pi - Pi,i8pi + VjS<f>i + [Pl, Xl + P^XrXj^j) (plHl+ 



+ (^H 2 + ^[Pj, Xj + Pj, XjX 



JPJ 



(3.13) 



8 Note that it suffices to consider a first order expansion in the energy and momentum constraints since 
higher order contributions vanish on substitution into the action. 
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The above can be substituted back into the action (3.5) expanded to the desired order and, 
after removing total derivatives and using the background equations of motion, yields the 
perturbed action in terms of the n physical degree of freedom 5(f>i. For simplicity we begin 
with the second order action 



P Xl + Px lXl cpj) 5$ - PxMiMi 



2Pl,IX I 4 , I$IJ ~ -TrPl,XiX I Pj,Xj4 > I^J 



no<pj 



(3.14) 



where the effective mass matrix is given by: 



Mu = -Pi,iiS u + v JJ + 7^ (Pi,ix i Pj > Xj4>Uj + Pj,jx j Pi,x i 4>'j4>i 

1 . ... l d /q3 \ 



(3.15) 



which can be used to evaluate the two point function of field fluctuations, as in (2.12), using 
the standard prescription. We then find, after some lengthy calculations, the corresponding 
third order action (see [49, 50, 62] for analogous calculations) 



5(3) = dtd 3 xa 3 



3H 2 a 2 + 2Hap]. + - [P^p] ■ - %/3 



a 



I\i I Pl,X! 



( fit® 2 ~ ^iaS4>! + ^5$ 



fafi + ±6<j> l f\ 5^ P ItXlXl Xi + lPi,ix lXl XiS<i>i 



( ^/« 2 - <j>iat8<pi + ^5$ 



^ + -8$ 



) Pi,ix z S<f>i + ^PijiXzXrStf 



+ -P IiIia 6<fi + -Pi,ni5<$ - Y WjjaSfaSfa - ^ £ ^V,ij K S<l>iS<f>j6</>K 



+ qPi,x i x i x i X i 



J K 



(3.16) 
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where for notational convenience we have introduced Xj = <frj y5(j)i — (j)iaj and we note again 

that only the first order energy (3.12) and momentum (3.13) constraints are required. The 
above results are consistent with an analogous calculation by [50], who consider a slightly 
more general action where P in Eq. (3.15) is a function of the kinetic functions Xjj = 
— \g^ v ^i^j^ and the scalar fields <f> = {0i, fa, </>n}- Note that where we use a~ 2 S^ to 
raise spatial indices, [50] use cP J only, which accounts for the additional factors of a~ 2 and 
a -4 in the latter's results. The two sets of expressions are identical when the metric is written 
explicitly. 

We now consider the leading contributions to the action (3.16) in slow variation and 
small sound speeds, since the dominant contribution in the following path integrals will be 
around horizon exit. Neglecting the purely gravitational part of the action and following the 
arguments regarding the more general versions of slow variation in [49, 50] we find 



'(3) 



dtd 3 x a 3 V] 



2 < PiPi,x I x I + -^4>I x iPi ,Xj,Xj,Xi 



\4>iPi,x I x I 54>i5cj) l I 54> I \ i 



(3.17) 



Consider now the multiple-DBI scenario, described by the Lagrangian (2.2), which we 
repeat for reference 

Pi = -j^(l- y/l - 2fWX?j . (3.18) 
It is then straight forward to compute the following derivatives 

Pi,x t = — ; -P/.X/X/ = P^XjXuXj = -i^-- (3.19) 

ci cj cj 

Upon substitution into (3.17) and keeping only terms at leading order small sound speeds, 
we arrive at 



5 { 3) ~ / dtd 3 x a 3 ^ 



1 1 

2^f 



1 1 

2^f' 



(3.20) 



This is the third-order action for field fluctuations in the multi-DBI scenario to leading order 
in slow variation and small sound speeds, which is justified around horizon exit. 
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3.2 The path integral formalism 



With the action (3.20) we are now in a position to calculate the three point function of field 
fluctuations at horizon crossing to leading order in slow variation and small sound speeds. To 
this end we adopt the path integral technique and, for brevity, we refer to [62] for a clear and 
detailed description of this method. We first require some standard results, the first being 
the propagator and its time derivative. The following is easily obtained from the second 
order theory (3.15) assuming slow variation at horizon exit, in exactly the same way as the 
power spectrum (2.12) 



H 2 

<^/ki(n)^jk a (T 2 )) = (2^) 3 ^g(l +tc/fciri) x 

(1 -ic / A;iT 2 )e- lfclC ^ ri - T2 ) ( 5 / > 3 )(ki +k 2 + k 3 ), (3.21) 

7a(l + iahT^e-^^-^SuS^i^ + k 2 + k 3 ). (3.22) 



where r represents conformal time. Indeed, by considering equal times t\ = t 2 in the super- 
horizon limit |c/&it| <C 1, the above yields exactly the definition of the dimensionless power 
spectrum (2.12). In addition we require the following time integrals, which can be obtained 
by choosing the appropriate contour in the complex plane (analogous integrals appear in [48]) 



Kci 7_oo {Kci) 1 J_ oc [KaY 



Given the above we are now in a position to use the standard prescription, as in [62], to find 
the following contributions to {^(pi^d^j^d^Kiia) from the first and second terms in (3.20) 
respectively 



^\^km^ SljBl ^ + k2 + k3) ' (3 - 24) 



HA i fcg(k 2 -k 3 ) /i (k 2 + k 3 ) 2k 2 k 3 \ (3) 

(27r) 471^7^ iPI \k + ^ 2 ~ + ^) 5ij6ik 5 { kl+ k2 +ks) ' (3 - 25) 
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where terms in the above are to be evaluated at horizon crossing, since this is the dominant 
contribution to the relevant time integrals. Note that there is a sign ambiguity in the above 
from using (2.6) to write Hj in terms of e^. The above are valid provided we assume <pi < 0, 
which will be the case in the scenarios we consider. Finally then, we sum these contributions 
to arrive at the three point function of field fluctuations at horizon exit, to leading order in 
slow variation and small sound speeds 



where the fc-dependent parameter A is given by 



A(k 1 ,k 2 ,k 3 ) 



^(k 2 -k 3 )(- + 



1 , (k 2 + k 3 ) , 2k 2 k 3 



K 2 



+ 



K 3 



+ perms 



(3.27) 



It is then trivial to check that the above expressions recover the single field result [79]. 



3.3 The equilateral non-linearity parameter 

Given the expression for the three point function of the field fluctuations at horizon exit 
(3.26), it is then straight forward to find the corresponding contribution to the three point 
function of £ using the 5N formalism. To this end, we substitute the result (3.26) into the 
first term in (2.13) and consider the two field scenario with fields <f> and \ 



= (*0'ffi A "fry ( ^= + -£L=) *0>(k 1 + k. + k.), (3.28) 

where the subscript 'eq' denotes that we are considering only the equilateral contribution. 
Comparison of the above with the definition of the bispectrum (2.8) gives, 



B ( (k 1 ,k 2 ,k 3 ) eci = - 



1 H* A(ki,k 2 , k 3 ) 



4 ci 



+ 



2ei x) c, 



4tt 4 P, 



J2 i kf 



1 A(fci,fc 2 , 




(3.29) 
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where we have used (2.11) to replace V* with V^. Finally then, the terms in parenthesis can 

f (3) 
'NX 



(3) 

be associated with the non-linearity parameter f^J by inspection of (2.15) 




A3),, , , x 5 1 A(fci,fc 2 ,fc 3 ) , 



The above expression for is the main result of this section and, in the absence of addi- 
tional dynamical restrictions to evaluate the derivatives of N, cannot be developed further 
analytically. The ^-dependence is unchanged compared to that of the single field scenario 
[79], such that this contribution does indeed peak in the equilateral limit k\ = = k^. By 
considering x — > 0, such that N tX — >■ and — > and using the background equations 

of motion (2.3), it can be shown that the final terms becomes equal to one, recovering the 
single field result. 

We notice then that the above has the form of the single field result, which is precluded 
observationally by the strong c~ 2 dependence, modulated by an expression dependent on the 
background evolution after horizon exit. It seems possible then that, in some circumstances, 

(3) 

this modulation may suppress the value of /^l to remain within observational bounds. Such 
modulation has been found in similar multiple-field DBI scenarios. For example, [49] use the 
adiabatic-entropy perturbation basis to obtain the single field result modulated by a cos 2 
term, where depends on the background trajectory after horizon exit. This modulation has 
then been exploited in concrete examples to suppress the otherwise observationally precluded 

(3) 

value of /j^jl, as in [60]. It is not a priori obvious if this possible in our case by inspection 
of (3.30) however, since the behaviour of this term is highly model dependent. We address 
this issue in the following section by considering scenarios in which the derivatives of N can 
be evaluated analytically, as we did for the local contribution in [59]. 



4 The total three point correlation function 

The expressions for f^) (3.30) and f$ (2.16) contain field derivatives of the number of e- 
folds N. Given the lack of a unique attractor in multiple-field scenarios however, additional 
restrictions to the background dynamics are required to further develop such expressions 
analytically. For example, in our previous work [59] we adopted the method of [29] and 
demanded a sum separable Hubble parameter 9 . Not only did this allow the violation of slow 
variation after horizon exit, it also suited the case of non-standard kinetic terms, since the 
dynamics are better described in the Hamilton-Jacobi formalism. Given this restriction, we 
exploited the resultant integral of motion to derive analytic expressions for the derivatives of 
N and in turn fj^ in the multiple-DBI case. In this section we apply those results to present 

(3) 

an analogous expression for which, to the best of our knowledge, is the first example of 
the application of the separable technique towards the equilateral contribution. 

9 As opposed to a sum separable potential V((f>,x) = V((f>) + V(x), as was used in the original application 
of this technique [21]. 
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To make analytical progress we now restrict our attention to two-field models, with 
fields (ft and %, that posses a sum separable Hubble parameter 



H( ( ft,x) = H^\<ft) + H^\x), (4.1) 

which leads to a number of simplifications. Inspection of (2.5) shows that the sound speed 
becomes a function of its respective field cfti only, such that c^\<ft) and c^ x \x)- Moreover, 
mixed derivatives of H (i.e. H^^) become zero. The combination of the above reduces the 
number of relevant slow variation parameters 



e 



W =2c W(i^M , e W = 2 C W ^) , (4.2) 



V W = 2c W%, r,^ = 2cW%^, (4.3) 

,W =2c W^ > s (x) =2c (x& (4 . 4) 



where e = e^) + e^ x ^ and we emphasise again that rp 1 ^ and s^ 1 ' can become much greater than 
one after horizon exit. Crucially though, the above assumption enables us to calculate the 
field derivatives of N analytically by defining an integral of motion. Here we simply quote 
the results of our previous work [59], where the full details of the calculation can be found. 
The derivatives of N can be expressed in terms of slow variation parameters 

= r== u, N, Xt = r^— v, (4.5) 
where, for brevity, we have introduced the following definitions 

u= * f , v= * f , Z f = ^^ !—f-. (4.6) 

U* ±1* £/ 

On substitution of the above into the result (3.30) , we arrive at the following expression for 

f (3) 



(3) 



5A(k 1 ,k 2 ,k 3 ) 1 



+ 



6 E« k f cl 



+ 



,(\) 



(4.7) 
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which we emphasise is valid for the two-field DBI scenario assuming comparable small sound 
speeds and slow variation at horizon exit, in addition to a separable Hubble parameter. By 
setting x ~ > 0, we find Zf — > u — > 1 and v — >■ such that (4.7) again recovers the single 

field result [79]. The above is analogous to the result (3.30), in that we find the single field 
result modulated by a term dependent on the background evolution after horizon exit. The 
advantage here however, is that the behaviour of the modulation becomes more transparent. 
For example, we note that this term is approximately 0(1) in slow variation at horizon exit. 
Moreover, since neither of the terms in the numerator are positive definite it is conceivable 
that, with a sufficient level of cancellation, the modulation term may suppress the otherwise 
prohibitively large contribution of c~ 2 , providing an observationally viable value of f& at 
the end of inflation. It remains to be seen if this is possible in a concrete setup however, 
which we intend to address in the following section. 

Alongside the above equilateral contribution (4.7), we must also consider the local signal. 
Here we can directly quote the analogous expression for fj^£ from [59] 



'NL 



in 



X4>) 



+ 



Xx) 



(r, 



(x) , „(x) 



Xx) 



v +2 



Xx) 



A 



( U 2 I V 2 
~W> t ~Xx) 



(4.8) 



where the parameter A is defined as 



and we have introduced 



(4.10) 



These additional parameters appear because the expression for (2-16) contains second 
derivatives of N, whilst (3.30) has only first. Note that, whilst the terms preceding the 
parenthesis in the expression for A (4.9) are 0(e*), r]j s and s s j s can become much larger than 
unity, producing observable local type non-Gaussianity. We demonstrated this in a concrete 
model in [59] by considering inflation in the tip regions of two warped throats, in which s^ 5 
becomes enhanced by the abrupt change of and at the end of inflation. 

Equations (4.7) and (4.8) together provide the full expression for /tvx, and in turn the 
bispectrum (2.15), to leading order in slow variation and small, comparable sound speeds at 
horizon exit, given a sum separable Hubble parameter. We have noted that both contributions 
can, in principle at least, produce contributions to the bispectrum, rendering this one of the 
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few explicit models capable of doing so (see [60, 63] for alternatives). It remains to be 
seen if this is possible in practice however. Ideally, a consistency relation between the two 
contributions would elucidate this point further but, given the considerable freedom within 
the model, finding a general relation has so far proved difficult. In the absence of a consistency 
relation therefore, it is useful to consider some specific models, as we did in [59] to study 
the local contribution produced by inflation in the tip regions of the throats. We leave this 
for the following section and here simply highlight that this scenario does indeed provide the 
potential mechanism to produce a mixed non-Gaussian signal. 

Before proceeding to specific scenarios, we first note that whilst the expressions (4.7) 
and (4.8) describe the production and evolution of non-Gaussianity during inflation, these are 
not necessarily the final observed values. For completeness we must track their evolution from 
the end of inflation until they are imprinted on the cosmic microwave background (CMB) at 
decoupling. Given our lack of knowledge of the early Universe however, this is not generally 
feasible. As such, recent work has considered whether such non-Gaussianity produced during 
inflation can indeed imprint upon the CMB [33-40]. Such study provides valuable clues as to 
what we can infer about inflationary dynamics from observations of non-Gaussianity and it 
would be interesting to include such considerations in this scenario. Here however, we simply 
illustrate the potential production of mixed non-Gaussianity through multiple-field dynamics 
and small sound speeds during inflation. 

5 Illustrative Examples 

In the previous section we demonstrated the possibility of producing mixed local and eqUi- 
ZoX 

lateral non-Gaussianity in the multiple-DBI scenario, through the expressions for f^J (4.7) 

and (4.8). The aim of this section is to look more closely at the feasibility of producing 
such a signal by considering some concrete examples. We begin by revisiting the case of 
inflation in the tip regions of two warped throats, where previously we found enhanced local 
non-Gaussianity caused by the rapid increase in the sound speeds at the end of inflation. We 
find that the required modulation is insufficient to suppress the prohibitively large contribu- 
tion of 2 in (4.7) however, and leave further exploration of parameter space to future work. 
We then consider a phenomenologically similar but fully analytic model by choosing expo- 
nential warp factors, finding a regime in which observationally viable mixed non-Gaussianity 
is produced. Whilst not derived from string theory, this does provide a first step towards 
understanding the regimes in which this occurs, before progressing to more realistic setups. 

5.1 Inflation in two cut-off throats 

Previously, we studied the case of inflation in the tip regions of two warped throats and found 
enhanced local non-Gaussianity at the end of inflation, caused by a sudden increase in the 
sound speeds [59]. Here we revisit exactly this model and include the equilateral contribution, 
to study whether the background evolution is sufficient to suppress the contribution of c~ 2 
in (4.7). We keep our discussion of the model brief, since full details can be found in [59]. 
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We choose a model whereby two probe D3 branes traverse two distinct warped throats 
glued to a compact Calabi-Yau in type IIB string theory [56-58]. As such, the warp factors 
are given by 

' w -M i+A2iog (£»' (5 - i) 

where is given by replacing — >■ x an d for simplicity we have assumed the same warping 
in each throat. A full discussion of the physical significance of the Aj and the infrared 
singularity at = A3 e -1 / A2 can be found in [59]. We also choose a linear, separable Hubble 
parameter 

H(^ X ) = U^U + U^x, (5.2) 

such that 0. Using (2.5), we arrive at the following expression for the sound 

speeds 



1 + ^ 1 + A 2 log U 



(5.3) 



where is again given by replacing <f> — )■ x- We then solve the background field equations 
(2.3) with Ai = 6 x 10 16 , which is typically required in standard DBI [47], A2 = 2 and 
A3 = 1. Furthermore, we choose </>(i*) = x(t*) = 1 as our initial conditions, where t* is 
the time at which observable scales exit the horizon. We require N — 60 e-folds between t+ 
and the end of inflation, which we choose to define as e = 1. Finally we introduce a small 
asymmetry such that VS® = 1.188 x 10~ 6 and = 1.192 x 10~ 6 . With this choice of 
parameters c^'/cW ~ 1 + 10 -3 at t* and is therefore consistent with our approximation that 
the sound horizons are comparable when observable scales exit. Given this, we solve the field 
equations and plot the inflationary trajectory in figure 1. Whilst approximately straight for 
most of inflation, the trajectory finally curves in the x direction as the sound speeds rapidly 
increase in the tip regions of the throats. The turn in the trajectory is due to our choice of 
%{x) such that c&) increases slightly before towards the end of inflation. 

Given this trajectory, we study the evolution of the relevant quantities as a function of tf for 
a fixed t*. We find that the quantities associated with the two-point statistics are consistent 
with observations [80], where = 2.44 x 10 -9 and — 1 = —0.0108 at the end of inflation 
(see [59] for the full expressions of these observables in the separable Hubble approach). 
With respect to the three-point function, figure 2 shows the evolution of an d /^l m 
the equilateral configuration, calculated using (2.16) and (3.30) respectively, as a function 
of the slow roll parameter e. As discussed in [59], the rapidly increasing sound speeds in 
the tip of the throats produce fj?) ~ —20 at the end of inflation. Looking now at the 
equilateral contribution, we see that the value of c* ~ 10 -2 provides /i? ~ 10 4 at horizon 
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Figure 1. Left: The trajectory in field space, originating at </>(£*) = x(i*) = 1 and ending after 
TV ~ 60 c- folds of inflation, when e = 1. Right: Enlarged region of the trajectory (solid) illustrating 
the turn in the x direction towards the end of inflation. The dashed line shows the straight line 
trajectory corresponding to = WW for comparison. 
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Figure 2. Left: Evolution of with respect to tf, plotted as a function of e for the trajectory in 

(3) 

figure f . The background evolution after horizon exit enhances the initial value, producing /^l ~ 
f2700 at the end of inflation. Right: The evolution of where the rapidly varying sound speeds 
produce /i^ ~ —20 at the end of inflation. 



exit. Thereafter however, the background evolution actually enhances this value between 

(3) 

horizon exit and the end of inflation, eventually giving f£{ ~ 12700. We conclude then that 
this scenario is inconsistent with current observations [80]. Whilst this need not necessarily 
be the case for all parameter values, we leave a fuller exploration of the parameter space to 
future work. 
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5.2 Exponential warp factors 



In this section we consider a phenomenologically similar, analytic model that is able to 
produce both equilateral and local non-Gaussianity that satisfies current constraints. Let us 
consider the following separable Hubble parameter 



H = Hn 1 



.4 



.4 



2 2 

We make the following choice for the warp factors 



x e -p x 



(<f>) 



f 



fix) 



1 



e 2( 7 -l)a^ 



(5.4) 



(5.5) 



(5.6) 



All the parameters in the previous formulae are positive. This choice of warp factors is 
motivated by the fact that it will lead to an analytically solvable system. The sound speeds 
read 



,00 



B(f> e -(7-lW 



The equation of motion for the scalar (j) results 

= -2c^H i<t> = —H a A<f, B^ e" 1 
where the solution, assuming (f)(0) = (/)*, is 



<f>(t) 



1 . 
H In 

a 7 



1 - ( c? Hq 7 Aa, Ba, e 



-7 a (pi, 



(5.7) 
(5.8) 

(5.9) 
(5.10) 



while an analogous solution can be found for x{t)- The solution for the scalar field is a 
decreasing function of time t. The speeds of sound are increasing functions of time, provided 
that 7 and 5 are larger than one. The slow variation parameters read (in the limit in which 

A 6 e~ a ^ < 1) 



^ = \^LB,o?Al 



M) 



2 H 2 
H 



-^B^a 2 A+e-r*, 
ti 



Ho 
H 



B^ ( 7 -l) a 2 ^e"^ 



and we note the useful relations 

M 



H 



-\/2 cW a \/eW 

-a/2 cW (7- 1) 
Hp 

.. _ ryW rjix) 
a 2 c (jt>) p2 c ( x ) 



a 



(5.11) 

(5.12) 

(5.13) 
(5.14) 

(5.15) 
(5.16) 
(5.17) 
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The number of e-folds is then given by 




Ne ~ f~ ^ I ™ ™ I f 2/32 (5-1) ( C (X) C (X) J (5 ' 18) 

(5.19) 



We would like to find an inflationary trajectory for which the amplitude of the non- 
Gaussianity parameter / N ^ can be tuned sufficiently small to satisfy present constraints. At 

the same time, we would like that the amplitude of the local non-Gaussianity parameter / N ^ 

is sufficiently large to be detectable in the future, but still satisfying present-day constraints. 

(3) 

Recall the expression for (4.7) where 



u = - + — e W - 6 W + y J_ - 1_ , (5.20) 
2 2e ^ a V^) p^x) )' 

v = 1-u, (5.21) 

(3) 

in the approximation e* <C 1. If the speed of sound c* <C 1, the amplitude of is pro- 
hibitively large, unless u is tuned in such a way that, at the end of inflation, the following 
inequality is satisfied 



„3 



+ 



(l-u) 



a{u) ee )* * ' « 1. (5.22) 

/2 (1 — u) 1 



Recall also the expression for / N ^ (4.8), where the dominant contribution is given by 

2 

2 ' 



where 



/nl = / I Tt4~A (5.23) 



( It 2 I V 2 



A= m ± rr\i + ^)- (5 - 24) 

In order to obtain a detectable /nl, we have to find situations in which either or both ?f s 
and s ss are large at the end of inflation. 

(3) 

Let us start discussing the conditions for being able to tune the value of / NL \ We assume 
the following hierarchy between the slow-roll parameters at horizon exit 
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ei 0) = rel x) , with r > 1. (5.25) 

2 2 

The function u{u), defined in (5.22), vanishes at the point r3/(rs — 1). We then demand 
that the parameter u satisfies 

2 

r 3 

Uend = ~2 \~ \ (5-26) 

r3 — 1 

at the end of inflation, where A small in absolute value. Expanding the function cr(u) at first 
order in A, one finds 

a(u end ) = -± — A. (5.27) 

We can therefore set cr(u en d) to be sufficiently small to compensate the enhancement associ- 
ated with the speed of sound in (4.7). 

Let us assume that, at the end of inflation, is much larger than 



e (x) = Kie (<t>) ; with ki « 1. (5.28) 
Moreover, we write the following expressions for the speeds of sound at the end of inflation 

C W = 4e W > c« = JeW, (5.29) 
or p z 

such that we can write 



Uend = 1 " J^— + ~ — ) (5.30) 

1 + Kl V K 2 K 3 / 1 + Ki 

Tuning properly the Ki then, this quantity can assume the desired value of (5.26). 

We now consider possible inflationary trajectories, with the specific requirements listed 
above, that lead to observationally viable non-Gaussianities of both local and equilateral type. 
We would like to determine which conditions we have to impose to the model parameters 
in order to satisfy all our requirements. For simplicity, we assume that the 7 and 5 are 
very large, so that we can safely neglect corrections weighted by inverse powers of these 
parameters. We proceed by discussing one by one the conditions that fix our parameters. As 
discussed in the previous sections, in all our analysis we make the hypothesis that, at horizon 
exit, = ci x ^ = c*. This implies 



(7 — 1) a a (7 — 1) \B 
The quantity (j)* is associated with the speed of sound c* at horizon exit by 
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In ) . (5.32) 



a (7 — 1) \ c 
The condition (5.25) gives 

= y/f^A x , (5.33) 

where we neglect corrections that scale as I/7 and 1/5 since, as we assumed above, these are 
negligible. Condition (5.29) gives 



/2 1 

Ab = —7 F ^7, (5-34) 

A X = —7 

Hence, combining with (5.33), we find 

§-*S- < 5 - 36 > 

In the approximation .ff ~ i?o, the slow-roll parameters at horizon exit are given by 



(5.37) 

c*, (5.38) 
c*. (5.39) 
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^3 
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/3 2 5 




fc 2 


c*, 







The dominant contributions to the number of e-folds depends on the terms evaluated at 
horizon exit and with the aid of the previous equations we find 

<5 ' 4o) 

The dominant contributions to the non-Gaussian parameters in this scenario read as follows 
(we write only their amplitude, and not the scale dependence) 

f (3)_ 5 1 ( r l/3 + l)( r l/3-l) 3 



(4) _ 2 (r 2 / 3 - 1) Kl (5f3 + nyya) 
/nl " rV3 7 "7 . ,/o-x2- ^ 42 J 



QK2 /3K3 
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Without providing an exhaustive analysis, let us consider a concrete set-up in which we assign 
the following numerical values to the parameters «j 



ki = 10" 3 , (5.43) 

k 2 = 10" 1 , (5.44) 

«3 = 10. (5.45) 

We demand also that 

4 0) ^ 4 X) * HT 2 , (5.46) 
to obtain a sufficient number of e-folds. With this choice, the value of u en d is 

Uend - 1 + 6 X 1CT 3 , 

implying that (see (5.26)) 

r ~ 2.3 x 10 3 + 5.8 x 10 5 A. 

Hence, assuming that A < 10 -3 (we will see that this assumption is satisfied in our set-up) 
(5.36) gives 

/3~2a. (5.47) 
The condition (5.46) then provides the relations 

c * * (5 - 48) 

implying that 5 ~ IO7. Let us consider the local non-Gaussian parameter f^: using (5.42) 
we find 

/nl - -3 x W' 5 a 3 5. (5.49) 

Since our parameters a and 5 are positive, this quantity is always negative. Choosing for 
definiteness 

10 5 

5 ~ —3 (5.50) 
3 or 

we get 

l/Sl = O(l). 

Substituting the information of (5.50) into (5.48), we get 

c* = 3a x 10~ 7 . (5.51) 

We are not allowed to choose too large values for a, since condition (5.50) would lead to 
small values of S, against our working hypothesis. We set a = 3, that gives 

c* = 10- 6 . 
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Hence, f*J is easily obtained from (5.41) by substituting our values of the parameters: 

fSl ~ -1.4AX 10 8 . (5.52) 

By tuning appropriately the parameter A to a value of order 10 -7 one finds 

\f§l\ * O(l0). (5.53) 

and its sign depends on the sign of A. Hence, this semi-quantitative analysis shows that 
within this concrete model we can obtain relatively large non-Gaussianities both of local and 
equilateral type. It would be interesting to perform a more complete analysis of this set-up 
to determine more precisely its predictions. Moreover, since the exponential warp factors 
were introduced for analytical ease, it would also be prudent to use these results as a starting 
point for further parameter exploration of the previous example (section 5.1), which whilst 
more realistic is less analytically tractable. 



6 Conclusions 



We analysed a setup of multiple-field DBI inflation leading to mixed form of primordial 
non-Gaussianity, including equilateral and local bispectrum shapes. Previously, we stud- 
ied a multiple-DBI model as an example of multi-component inflation with non-standard 
kinetic terms, showing that rapidly varying sound speeds can produce large local type non- 
Gaussianity during a turn in the trajectory [59]. Here we have included the equilateral 
contribution produced on sub-horizon scales, and found the possibility of an observationally 
viable mixed non- Gaussian signal. 

We used a general formalism based on the Hamilton-Jacobi approach, allowing us to 
go beyond slow-roll, combining the three-point function for the fields at Hubble-exit with 
the non-linear evolution of super-Hubble scales calculated using the SN formalism. We were 
able to obtain analytic results by taking a separable Ansatz for the Hubble rate. We found 
general expressions for both the equilateral and local type non-Gaussianity parameter /nl- 
The equilateral non-Gaussianity includes the usual enhancement for small sound speeds, 
but multiplied by an analytic factor which can lead to a suppression of this quantity. This 
modulation may suppress the value of /^l to within observational bounds (see [49, 60] for 
analogous results). 

We applied our findings in two explicit scenarios. In the first model, previously found to 
have detectable local non-Gaussianity, we found that the equilateral signal is not sufficiently 
suppressed to evade current observational bounds. In our second set-up we constructed an 
observationally viable model which exhibits both equilateral and a negative local f^, 
providing a first step towards understanding regimes in which this can occur in more realistic 
scenarios. 

Open issues remain however, the most prevalent being a better understanding of the 
likelihood of such a mixed non-Gaussian signal. For example, it would be interesting to 
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use the results of our analytical example in section 5.2 to more methodically assess the 
parameter space of the phenomenologically similar but more realistic scenario in section 5.1. 
Moreover, a consistency relation between the more general formulae for (4.7) and 
(4.8) would be desirable, since this would be more widely applicable. The use of a sum 
separable Hubble parameter, whilst providing analytically tractable expressions, also restricts 
the choice of potential. Thus it would be more general still to consider alternative analytical 
or numerical methods. Finally, we again note that our approach only treats multiple-field 
dynamics during inflation and we cannot necessarily conclude that such values are observed in 
the CMB. It would be prudent therefore to use analogous techniques to [33-40] to further our 
understanding of inflationary dynamics from observations of non-Gaussianity in the CMB. 

We conclude by noting two further directions to develop the above, in addition to 
increased generality. In our previous work we made a first attempt at addressing the next 
order statistic, the trispectrum, by plotting the evolution of the analogues of the non-linearity 
parameter, <7nl and tnl ; in the example of inflation in two cut-off throats. In light of 
impending observations by Planck [5] however, it would be interesting to further characterise 
the model by more rigorously considering its predictions for the trispectrum. Finally, we 
note again that in the third order action (3.20) we worked to leading order in small sound 
speeds at horizon exit. It is conceivable however that the next to leading order terms may 
still provide a significant contribution, with the potential to generate shapes distinct from 
the local and equilateral types. Since tentative signals for the orthogonal shape have been 
detected by WMAP [81], and again given the impending results from Planck, it would be 
interesting to consider such contributions in future work. 
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